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Inleiding 


Het probleem van da vorming van rustende peotubevenwen in de corona, als 
gevolg van het optreden van een thermische instabiliteit, laat zich onder= 


verdelen in een aantal deelproblemen: 
1) Onder welke omstandigheden in de corona treedt een neiertache instabili= 


teit op? ‘ 
2) Hoe veranderen dichtheid, edt uitstraling en druk gedurende het 


condensatieproces? 
3) Wat ís de eindtoestand? 


Hier zullen we ons voornamelijk bezighouden met vraag 2), hoewel ook vraag 


1) aan de orde zal komen. 


Vraag 2) zou ook het antwoord op vraag 3) moeten leveren, Een aantal aanna= 
men echter (zoals: atmosfeer optisch dun, materie volledig geïoniseerd) ie 


niet meer juist voor de eiudtane. van de ontwikkeling van de protuberans. 


Er zijn berekeningen gemaakt voor condenserende gebieden van geringe afme= 
ting (lD=berekening, + 250 km) en voor grotere gebieden (2D-berekening, 

+ (15000 Xx 6000 km)). De afmeting van het condenserende gebied en de begin= 
dichtheid bepalen of de condensatie isobarisch verloopt. Bij grote afmetin=. Ne 
gen van het koelende gebied kan een aanzienlijke drukdaling in het centrum 


gedurende een groot deel van het condensatieproces optreden, hetgeen ook he 


kan plaatsvinden bij kleinere afmetingen, doch den bij een hogere begindichtheid. 


Bij de berekeningen is de invloed van de gravitatie nooit meegenomen. De 
rol van het magneetveld in de gemaakte berekeningen 18 vrijwel BAERERE tot | 
het effekt op de warmtegeleiding. | he 


Voor de berekeningen werd gebruik gemaakt van een vast rooster en het nume= 


rieke schema SHASTA (zie appendix en 2434). 


Soortgelijke berekeningen werden uítgevoerd door Hiláner (1974) . Hier wor= 
den echter enkele nieuwe punten aan de orde gebracht, zoals het piach= 
effekt (zie blz, !4 ), de invloed van de gebruikte stralingsfunktie (Cox & 
fPucker of McWhinter) en de wartragende invloed die de warmtegeleiding op de 
_ thermische instabiliteit kan hebben. Andere verschijnselen, zoals line a 
tying (zie lijn 12) en de tearing mode instabilities in een neutrale laag 
(zie 8), evenals de stabiliteit van de protuberans in het gravitatieveld Re 


van de zon (zie 9) worden hier niet beschouwd. 
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De thermische instabiliteit in de corona 


De corona bestaat uit een heet, ijl plasma, Er zijn drie mechanismen die 


het plasma doen stralen: 


i) vrij=vrij emissie (Bremsstrahlung); 

2) lijnemissie; 

3) vrij-gebonden emissie. 

In het geval van kosmische abondanties (zie 1) ziet de uitstraling als 
functie van de temperatuur er uit als in fig. 1 (C & T). Hier is uitgezet 
E(T) met f,(Te) Pe = uitstraling per seconde per volume eenheid. Merk 
op dat de uitstraling evenredig is met de (dichtheid)? (zie D.P. Cox & 
W.H. Tucker, 1969). 


Bij het berekenen van deze curve is een aantal veronderstellingen gemaakt: 


a} het gas is cptisch dun voor alle straling (d.w.z. foto-ionisatie ís te 


verwaarlozen en alle uitgezonden straling ontsnapt). 


b} alle geëxciteerde ionen deexciteren naar de grondtoestand in een tijd 


kort vergeleken met de tijd tussen twee botsingen, 
c) de karakteristieke tijden en afstanden waarover de dichtheid en de tem= 
peratuur aanzienlijk veranderen zijn lang vergeleken met de tijden die 
nodig zijn om een stationaire toestand te bereiken. 
De curve vertoont een maximum bij 2 . 10° K , vooral veroorzaakt door 
meervoudig geïoniseerde zuurstof. Een maximum bij 2. 10* K wordt ver- 
oorzaakt door Lya . Vrij-vrij emissie is alleen belangrijk voor tempera- 
turen boven de 107 K . Merk op dat het verloop van de stralingscurve (met 
aane de hoogte en de posítie van het maximum) zeer sterk afhankelijk is van 


de gekozen zuurstofabondantie. 


Daar het gas optisch dun is verondersteld, ontsnapt alle straling. Dit gaat 
ten koste van de inwendige energie van het gas, m.a.w. het coronale gas 

(temp. + 10° K) koelt snel af (zie de steïle tak van de stralingscurve bij 
10° K). Dientengevolge neent de schaälhoogte van de corona ook af, totdat 
het gas optisch dik wordt en er evenwicht ontstaat tussen uitstraling en 

absorptie. De corona is dan verdwenen. Dat de corona in werkelijkheid een. 
permanent verschijnsel is, heeft meerdere corzaken, zoals de dissipatie van 
nechanische energie, warmtegeleiding en de aanwezigheid van een ingevroren 


magneetveld. Deze drie punten zullen in het kort hier worden behandeld. 


De toevoer van mechanische energie is een gevolg van de dissipatie van gol- 
ven, Hierbij wordt de bewegingsenergie van de deeltjes die aan de golfbewe- 
ging meedoen, omgezet in inwendige energie (d.w.z. warmte) . 

De golf zal beter dissiperen naarmate er meer deeltjes per volume eenheid 


aanwezig zijn in het medium waar de golf doorloopt. De eenvoudigste veron” 


derstelling is dus, dat de mechanische energiedissipatie per volume eenheid 
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evenredig is met de dichtheid ps: 


rect en 8 
Hierin is C een functie die een weinig van de temperatuur T zal afhan= 
gen. Als T hoger wordt, d.w.z. als de random bewegingen van de deeltjes 
toenemen, zal de weerstand die de golf ondervindt, en daarmee de eea all 
tie, toenemen. Hier zuilen we C constant veronderstellen. 


De corona is nu energetisch in evenwicht wa we C zo kiezen dat: 
— £(T) p? +Cp = 0 


C = Ê(T) Pp. 
De corona is nu echter nog niet stabiel voor storingen in de dichtheid. Dit 


volgt uit de energievergelijking, door toepassing van een lineaire stabili- 


teitsanalyse: 
c Em - ET) p? + ECT) De P 
rot e “e 
Ed inf De temperatuur T iin 
Bierin zijn { De de { dichtheid  Y de corona in evenwicht en { p zi 


respectievelijk de hehehe in het gebied met de storing in de dicht- 


heid. Cp is de soortelijke warmte per volume eenheid bij constante druk: 


C, “ CoP (en per gram). Tevens geïdt voor een isobarische storing vol- 
gens de ideale en dat PT = Pe Tg waarmee de energievergelijking 


zich jaat schrijven als 
3T | … £@) 
Cp De Pa | Bar Te + ÊCT) 


T 
1 e eq 
met T = Te + T, 5 ir ak } , laat zich de energievergelijking na 


Taylorontwikkeling van £{(T) schrijven als 
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en r de storing in de temperatuur op t=0 is, 
> 
We zien dat de isobarische dichtheidsstoring instabiel wordt als a >0, 


d.w.z. 
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Áls in de buurt van T e= T. geldt É(T) se ze ‚ is deze conditie equiva= 


lent met 


| a < li | 


| 


Voor T= 10° K is in het geval van de Cox & Tucker stralingskromme 
ai , zodat de corona inderdaad instabiel is voor storingen in de 
dichtheid, 

Als de storing in de temperatuur op tijdstip OQ CT negatief is, zat 
Ere steeds negatiever worden. Dit komt overeen met afkoeling en met een 
voortdurende toename van de dichtheid, Daar waar Ti o > 0 vindt voort- 


durend verhitting plaats met een afname van de dichtheid. 


Dus als aïleen mechanische verhitting en straling de energiebalans van de 
corona zouden bepalen, zou deze spoedig opsplitsen in een aantal hete, 
ijle en een aantal koele, dichte gebieden en oock dan zou de corona ver= 


dwi jnen. 


Ex is echter een proces dat een isobarische storing in de dichtheid tegen= 


gaat, en wel de warmtegeleiding. 


Indien men de warmtegeleiding meeneemt in de energievergelijking, Îuidt 


deze: 
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waarin K de warmtegeleidingscoëfficient K = K el is. K is vrij= 
wel onafhankelijk van de dichtheid en de temperatuur. 
Laat L een karakteristieke afmeting van het gebied waarover de dicht- 


heidsstoring zich uitstrekt zijn, Dan kunnen we de warmtegeleidingsterm 


benaderen met 
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Lineariseren van de vergelijking levert: 


dT, [ ern mett 
en en = FIT en 
°p òt Pe Î 5 (TJ) 5 Te 
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De corona blijft nu stabiel als de afmetingen van het gebied met de dicht- 


heidsstoring voldoen aan: | 


5/2 
Od KE / 








{1 — Q) < 0 


Te 2 ie 
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of 


K, jr 7/2 b 
& 





L< L, met Lb, = 
(a) pe ET) 


Bij afmetingen groter dan L, is de warmtegeleiding niet meer in staat om 
de afkoeling (d.w.z, het inzetten van de thermische instabiliteit) te ver- 


hinderen, 


Bovenstaande conditie kan worden omschreven als 


Tmc STe 
met 
…Ì 
Í[ K, r 5/2 
Tae * | ad & karakteristieke tijd waarin de warmtege= 
Ne en Pa J leiding de storing T, Q teniet zou doen 
n ’ 
bij afwezigheid van koeling door straling. 
l Pe te 





Te = Ter ” mvm ola karakteristieke tijd waarin de 

: Ì T id | 

(tra) £( e) Pe storing T, , met een faktor e 
‚0 

zou aangroeien bij afwezigheid van 


warmtegeleiding. 


Tac 5 Te komt overeen met L= Le - _ 
In de corona, met Kk, ® 3.5 Xx 107S erg KTÌ sect emt (zonder magneet 
veld , of evenwijdig aan de veldlijnen). | OO 


T == 3106 « 
e 


an, © 0,35 # 10° deeltjes/cm° 


f (TT) Rs 10722 5 a == (Cox & Tucker functie; 


Omey * 0-6 McWhinter functie) 








geldt: L, Ss 78000 km. 


Gebieden in de corona waar de dichtheid wordt verstoord over een gebied 
ter grootte van ongeveer 80 000 x 80 000 km? re 2s 5200 sec tT) 
zijn dus stabiel, Gebieden met grotere afmetingen zijn thermisch insta- 
biel. 


Als gevolg van de aanwezigheid van magneetvelden in de corona, vordt een 


sterke anisotropie geïntroduceerd: Ki = 2, 1076 T5/2 


n? 
n in=l6 __© 





2 K) Á r°/2 n B … Ì Gauss 
rr = |E = 10°. | « 3 , 10° als Te = 10° K 
L, 1 he 


De + 0,35 « 10° / cma° 


De warmtegeleiding loodrecht op het veld kan dus worden verwaarloosd 
t.o.v. de stralingsafkoeling, mits de afmetingen van het beschouwde ge= 
bied voldoen aan: | 





S. 300 meter . 


Verdere informatie over dit onderwerp is o.a. te vinden bij E. Smith and 
E‚R. Priest (1977) en P,K, Raju (1968). 


Numerieke berekeningen 


De op te lossen vergelijkingen in vektorvorm vindt men in de appendix 

(zie blz. Al). Ze zijn hier în een vorm geschreven, die het mogelijk maakt 
het numerieke schema SHASTA er op los te laten (zie blz. Aá van de appen= 
dix). Alvorens we de resultaten van tweedimensionale berekeningen beschou= 


wen, kijken we eerst naar de uitkomsten van een aantal Eéndimensionale be= 


rekeningen. 


Eendimensionale berekeningen 


De snelheid van het geluid in de corona is ongeveer (bij een temperatuur 
van 10° K) gelijk aan Vv, = 1.1 x 10’ cm/sec. Deze snelheid is vooral 


van belang voor de drukverdeling. Zij  L de karakteristieke afmeting van 





aen gebied waarover een storing in de druk is aangebracht. Dan zal deze 

ë . 8 Ì 5 L ee , * 
storing teniet worden gedaan in een tijd aen. (Daarbij is aangenomen 
dat er geen magneetveld is en dat de drukstoring niet versterkt wordt door - 


andere oorzaken), In tabel 1 is deze tijd voor enkele afmetingen van L 




















berekend, Tevens is de karakteristieke koeltijd: Ta if 5 ” 
vermeld voor enkele dichtheden. | ‘ ak ECT) Pe 
Tabel 1 
L (km) Te (sec) De (deeltjes/cm°) Tg (sec) 

10 0,08 10° bh 160 

102 0.85 dû h 416 

10% 8,5 05° hh2 

10 85 ok} bh 

105 850 10 4.4 

105 8509 Kaal 0.44 

12 
En 3 Te * El (als Uu = gemiddeld molecuulgewicht = 1) 
8 Te | 
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De tijd Te uit bovenstaande tabel is niet meer dan een indicatie voor de 


afkoeltijd. Bij de berekening van T, is uitgegaan van een constante uit- 


C 
straling (eigenlijk: netto energieverlies per sec per cm°) RAD = É(T) Pp 


gedurende de afkseling. KAD zal veel te groot gekozen zijn als het plasma 
als gevolg van een zeer kleine dichtheidsperturbatie begint te stralen. 
Echter de laatste fase van de afkoeling kan zeer snel verlopen, daar RAD 
dan veel groter zal zijn dan ECT, = 10° K) De Met andere woorden de in 


de tabel opgegeven T, kan veel te klein zijn als we uitgaan van een zeer 


C 
kleine dichtheidsperturbatie, doch veel te groot in het geval van een 


dichtheidsperturbatie van de orde van grootte van de dichtheid zelf, 


Bovendien gaat de gebruikte formule uit van een isobarisch proces. In het 


geval dat de afkoeling zo snel verloopt, dat de druk zeer snel daalt, is 


| P Ed f @ Ld 
To Fogg en , Lmmers alleen de thermische energie van het af 


nd Y == À al 
f (To) Pe 


koelende gas is nog maar van belang, terwijl de energie van het instromende 


C 


gas kan worden verwaarloosd. T kan daardoor met maximaal een factor 
Y = z worden verkleind. De werkelijke situatie tijdens de snelle afkoel- 


Gaaf Cn 


fase zal hier tussenin liggen: het proces verliest zijn isobarische karak- 
ter, hoewel de compressie-energie die vrijkomt als het hete gas de conden- 


satie binnenstroomt, alsmede de afkoeling van dit hete gas, wel degelijk 


Panna oe Ee ann ne een teer 4 er ved mek rn ee Bike on ae Ass rade ne ee dee 
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van belang is voor de afkoelsnelheid van de coridensatie. 


Hoewel bij de hoge elektroaendichtheden (a, = 1052 8 101% cn”) het 
gas niet meer zal stralen volgens de Cox & Tucker functie, is het toch in= 
teressant om te zien hoe het proces verloopt als men toch de Cox & Tucker 
etraïiingsfunetie gebruikt, De afkoeling verloopt dan nameiijk zo snel dát 
nog voordat de materie rondom het afkoelende gebied sterk versneld is ou= 
der invloed van de drukgradiënt, de afkoeling reeds heeft plaatsgevonden 
in het gebied waar op t = 0 de dichtheid verstoord werd. We hebben hier 
dus te maken met een vrijwel isodensiech proces, De uitstraling als func= 
tie van de tijd volgt de Cox & Tucker functie van hoge naar Lage tempera 
tuur, Een voorbeeld hiervan is te zien in fig. 2 a t/m f, De dichtheid in 
het niet geperturbeerde gebied is Pe = 307Ìl g/em* , de perturbatie in de 
dichtheid (over 72.5 km) bedraagt 5% op tw 0 , Eet totale gebied dat bi: 
de berekening is betrokken heeft: een afmeting van 257,5 km . De mechant= 
sche verhitting is hier constant in de tijd verondersteld. De situatie op 


te Q is aangegeven in onderstaande figuur. 
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RE p ä én ” 8 a” en é 3 
Situatie op t e= 0 met drukevenwicht: pp, = 1 El g/en 


hond 


tT, = 10° K 


Kd 


De opgelegde randcondieties laten vrije instroming door de randen toe 





| 8 | 
* ® « s | P 3 
Yocr dit geval 15 Be 0,65 sec en To” te pl Seles 
(CT ) p 
. e Ts an 
R MECH RAD 3 d j 
op £=0 ie QE PP papel 1,052 - 1E} ETD 0 = 
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e 
Dat het proces niet isobarisch verloopt, blijkt wel uit de drukprofielen 


(Fig. Ze), en daar de dichtheid vrijwel constant is gedurende de afkoeling, 


“ie | : 


vak uit de tenperatuurprofielen (Eig. 2). Zen snelhetdspiek ontstaat daar 
waar het koele gebied overgaat in het hete gebied (fig. 2e). De afgeplatta 
vorm van de dichtheidspiek (fig. 2d) js een gevolg van het clipping effekt 
van SHASTA (zie appendix 


Het besproken geval wordt gekenmerkt door een zeer grote drukgradiënt op 
sen vrijwel vaste positie. Voert men cok een warmtegeleidingstern in, dan 
ordt deze drukgradient minder steil, terwijl bovendien de plaats waar deas 
Erukgradi ant optreedt, verschuift. Bovendien verloopt de «fkceling langzsa= 
mex. In fig, 3 a t/fme is het vereenvoudigde geval van een tempereatuurput 
(temperatuur in de put: 950 000 K, buiten de put: 10° «) geschetst, zonder 
stralingssfkoeling en met een warmtegeleidingscoëfficiënt K 6 in (re 
“2. 10° e‚g.s.). pam 1,67 % 10715 giem’ $ Er is deukevenwicht op te 0, 
Tijdens het proces is er geen drukevenwicht, daar de warutegeleiding zaer 
groot is, De karakteristieke tijd voor warmtegeleiding is ongeveer 
(L = 87,5 km} : 


T as 


HC DR , a 0,008 eec , 


Aaf bo 
si 


Dit ie dus ongeveer de tijd waarin de tewperatuurput wordt opgevuld, In de= 
ze tijd {0,008 sec} legt een geluidsgolf een afstand vaa } km af , waaruit 


blijkt dat er inderdaad gaen drukevenwicht kan zijn, 


a 


det gavolg van het warmtetransport is, dat de temseratuur in de tempera 
5Ci 


55 


tuurput stijgt, terwijl de temperstuur erbuiten daalt (fig. 3a). Hetzelfde 
gebeurt met de druk (Fig. 3E), zodat de materie versneld wordt in een 
richting die wan de put afwijst, Zelfs rudst da temvervratuurgrvadiënten hele 


$ 
5e 
maal verdwenen zijn {na ongeveer 0,02 sec} is de enelheid pas opgelopen tot 


3,3 X 10’ eu/sec , een Fractie 3 x 107* van ae geluidssnelheid, 


Î W, | ne Ì 
Laar cp £ = 0 ongeveer z van het beschouwde gebied heet en 5 koal 


was, bedraagt de eindtemperatuur ongeveer: 


} EO 
LT, 2 5. IOS +. 950 000 = 983 900 K 
£ 3 3 
Keren we nu terug naar het probleem van de stralingsafkoeling bij hoge 


dichtheid (p = 10711 g/cm®) met warmtegeleiding XK = 107* r5/2 (zie fig. 
6 a t/m e), expliciet meegenomen in d rgievergelijking (zie appendix) 
Ee), exp 2 genomen In de energievergelijking : APP bn 
dan zien we dat de afkoeling trager verloopt en dat în korte tijd het ge= 
hele beschouwde gebied afkoelt, immers de warmtegeleiding (hier met karak- 
teristieke tijd van ongeveer | sec} start het straliugsafkoelingsproces 
ook in die gebieden, die op t= G nog heet (10° K) zijn. De opttedende 


ge 
sneïiheden (fig. 4e) zijn veel lager dan in hef geval zonder warmtegelei- 


19. ier 


ding, terwijl de snelheidsprofielea ook veel breder zijn, hetgeen zijn ver 


klaring vindt ín de veel kleineze drukgradiënten over een groter gebied, 


Bij een veel Lagere dichtheid (5.10° eu”°} verloopt het proces gehesì an” 
ders, hetgeen vooral duidelijk tot uitdrcukking komt ia de kromme die hat 
nettovenergieverlies per voluwmeveenheid per rijdseenheid (RAD) HesEBenSE, 
(Éig. 5). | 
Op het eerste gezicht lijkt de kromme veel op die bij hoge dichtheid: eerst 
een sterke toename van RAD (XE), daarna een veel langzamere toename AKE), 
gevolgd door een zeer steile piek (III). In hèt geval van de hoge dichtheid 
herkennen we direkt de Cox & Tucker stralingskromme: 

L 650 000 K — 600 000 K 

IL 650 690 K — 325 GOO K 
KIT 325 000 K — 25 000 K ; piek bij 213 


1e 


oe 


Er 


en dr 


(Bij 20 COG X is nog eeu heel klein piekje van RAU nauwe: Lijks bespeur” 
baar.} 
Tu het geval van de lege dichtheid, dus als het proces veijwel isobarisch 


werlcopt (afgezien van de snelle fase, zie later}, wordt de invloed van de 
dichtheid up de uitstraling over een groot temperatuurbereik sterker dan de 


Suvloed van de stralingskromme die het temperatuurafhankelijke deel ven de 


uitstraling bepaalt: 

Ï 200 200 K — 60 000 K 

Iï 60 500 K — 25 600 K 
EEL 25 OO K — 10 GOO K . | | 
Dus hoewel aa het passeren van het maximum van de Cox & Tucker » kromme de 
functie £(T) sterk begiat af te vallen, i3 de toename van de dichtheid 
aanvankeïijijk meer dan woldoende om dit op te vangen, waardoor RAD neg 
blijft toenemen totdat een. temperatuur van ongeveer 60 200 K is bereikt. 
Daar vindt een kleine inzinking plaats, Echter de zeer kleine piek van de 
Cox A Tucker = kromme bij ongeveer 20 000 K {Ly a ‘) veroorzaakt de 
steile piek van RAD, De waarde van RAD bij 20 000 K is ougeveer òX z0 
groot als bij 200 GOO K , immers, de stralingsfunctie van Cox & Tucker 
dà 


met sen factor Ì coeneemt (zie fig. 5). 


hendtennotnnds af wet een factor 6, terwijl de dichtheid 


m3 
Cr 


valt GVEe 


Een tweede effekt dat optreedt is de pictselinge temperatuurstijging vlak 
nasst het koele gebied. Dit is eveneens in fig. 5a aangegeven: TEMP,, » 
Aanvankelijk, als de koeling nog langzaam verloopt, gaat de temperatuucda= 
ting gelijk op met die van het puat in het centrum van het koeïiende ge” 
bied, en wel omdat ook hier op t= 0 de dichtheid met 10% was verhoogd. 


Echter, als de afkoeiing in het centrum sneï gaat verlopen, begint de tem 
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jd . Deze verhouding schouwelt tussen 0,47 
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is vitgezet als functie van de ti 


ER ESM ET EE TEE ei ink innam: mnd nd nd hentai te. in naa me 
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ax & bd Ld ek  & 
en 1,45 ED is per definitie ĳî op t = 3á sec) , terwijl V zaif 
warieert van 3,51 « 10° am/sec tot 3,91 # 106 cm/sec en RAD van ! 
lop t = 0) tot 1710 . Kortom: Voaxtt) is vrijwel evenredig met 


RAD (Et) … 


Hat is ook interessant om enetheidsprofialen op verschiliende cijdetippan 
e bekijken (zie fig. 6, waar ook temperatuur” en dichtheidsprofielen zijn 
te zien). De materie begint te bewegen daar waar het gebied met de gepex” 
turbeerde dichtheid begint, immers hier ontstaat ean kleine drukgredient 
als de materie begint te koelen, Vervolgens begint cok de materie daarbuiT 
ken te bewegen met ongeveer dezelfde znelhaid. Na aen tijd gr EN 
à,5 sec (= de tijd waarin een geïuidsgolf van het centrum naär de rand — 
beweegt) begint ook de materie aan de rand te bewegen. Zolang het proce 
isobarisch verlcopt, en dit ís het geval totdat de temperatuur tof cage 


veer 200 000 K ie gedaald, is er nauwelijks verschil tussen de snelheid 


dichtbij het koelende gebied en aan de rand van bet rekengebied. Bij nog 
lagere temperaturen wordt de snelheid nabij het koelende gebiad aanmerka= 


ijk groter dan ean de rand, 


Beschouwen we vervoigens ket bocipsecek bij Aaneh enen hd een magneet” 
ctan COX, 


erts oe A % kad Q aan” 
En Ir i | 
respondeert met de neutrale laag, waar B= 0 5; B, *= 5 Gauss ; 


d pe sen e é | dak 
a = 1.44 * 1077 emt} . Op t=0 is de temperatuur overal 10° XK, 


veld, Het wagneetvelid heeft de vorm B » B 


L er drukevenwicht is, dat wil zeggen de totale druk (magnetische 


ss @ 


wij 
druk + gasdruk} ís voor alle punten gelijk. Zij 8 de verhcuding 


asdruk | 
gg = ti A n 
magnetische druk a? 


Br 


N Ed s Ld e “% k 4 @ dad kad $ 4 : REN 
De deeltjesdichtheid ver buiten de neutrale laag 18 1C” cu en is in Ùs 


utraie laser ongevee kee ho, zodat = g | | ek 

neutrale Laeg ongeveer 8 keer zo hoog, zodat B, Pr eutrale laag 

B = 0,14 3 B: = 543 3 B} op 175 km van het cenkrub. De gasdruk de 
Pi 


ireert dus over de magnetische druk bianen een afstand van ongeveer 


2) 


175 km van de neutrcaie laag. Dichtbij het centrum waar 8 > Î heef: het 


magneetveld aanvankelijk dus weinig invloed op het verloop van het koe! 


od Ed 


proces, Vanwege het hoge geieidingsvermogen. (hier is U = se veronder- 

steld} worêt echter steeds meer magnetische flux meegevoerd tot dichtbij 
de neutrale laag, waar de veldsterkte dan cok sterk toeneemt (zie fig. 7, 
Op t= 16,/ sec naren B even buiten de neutrale laag: 


Bo ‘t = HIG 1,48 (B, PS Na 11.62) . De invloed van het magneetveld 


ü 
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GN 
md 
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is dan ook van belang voor de dynamica van het gas: drukverschillen die 
ontstaan door afkoeling vond de neutrale laag, worden niet alleen maer ze 
goupensserd (het proces ja vrijwel isobariach}), door instroming van Zas, 
doch ook door een toename van de magnetischa veldsterkte. 
Ean tweede effekt van de toeneme van de veldsterkte zou kunnen zijn dat de 
tearing mode instabiliteit zou kunnen gaan optreden (zie Kuperus, M, aad 
haben E,, 1967). Hier gaan we daar echter niet verdar op in, 
Een derde verschil met het geval zonder magusstveld ia de vorm van de 
snelheidaverdeling kort na het begintijdetip t = © {zie fig. 7), Deze 
verloopt veel gladder, Dit is een gevolg van de eis van drukevenwicht op 
t = 0 bij een overal gelijke temperatuur van 19° K , die sanieiding tot 
glad verlopende dichtheid geeft. Daar per definitie het gebied aan de rand 


4 8 


geen netto energieverlies heeft, ja er op t = 0 en daarne een netto 
energieverlies in alle andere punten, en wel een groter eoergieverlies 
naarmate men dichter bij de neutrale leag (waar de dichtheid ket hoogst 

ie) komt, | 
Uiteindelijk treedt cok hier het pinch effekt op. Halveert men op t e= 89 gen 
de stapleugte DX, dan wordt de uiteindelijke afmeting van het koele gebied 


ook gehalveerd en is de en: wat hoger en vok de gradient van 5 





wordt groter (zie fig. 1Ì. 
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Tweedimensi onále bereken: „ngen 


bns tichendind ki enmecencieee Cane CROS EERDERE NO RATED GE EDGE: 


Een tweedimensionale geometrie biedt de mogelijkheid om anisotropie te ia= 


Ereduceren, Daze anisotropie ken bijvoorbeeld tot uiting komen ia 


„ de afmeringea van het afkoelende gebied; 
A E n , >} oe 
= het verleeop van de toestandsvariabelen (p , v , Bs, T , 5D) « 


In twee dimensies is het bovendien mogelijk om de zwaartekracht in te voer 
ren {verticaal} met zijdelingse aanvoer van materie edere Dit pros 


bleem worde hier echter niet (aumeriek} behandeld, 


Her plachreffekt 
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Het afnemen van de afmetingen van het afkoelende gebied vindt ook in twes 
dimensies piaat a, We zuilen dit verschijnsel hier het pinchreffekt noemen. 


Ld 


is o,a, te zien aan de dichtheidscontouren op verschillende tijdstip” 


ri 


UN en 4o 
Ge Le 


ven, Als voorbeeld beachouwen we fig. Ba, hb, oc, die de dichtheidscontouren 


weergeelt voor een afkoelend gebied dat op t = 0 de afmetingen 6DX X SDZ 
2h OGG X 5000 ke heeft. Uiteindelijk worden de afmetingen van het koelende 


gebied 3DK Xx 3DZ 12 000 x 25CO km , hetgeen wordt bepaald door de groot 
te van de stapleagte in de beide richtingen (reep, zerichting = verticaal: 
DX , ea zerichting = horizontaal: DZ). Pe dichtheid in dit gebied (3 x 3 


« 


punten) îs vrijwel uniform ala gevolg van het clipping effekt, dat ean ai= 
tenechap 1e van het gebruikte numerieke schems SRASTÁ, Door het opgeven van 
de stapgrootten DX en DA legt men dus de afmetiogen van bet koele gen 


a 


bied ua het optredes van de pinch, zelf vast, (Zie ondevataande figuur). 


DX Fia 
| Afmetingen uiteindelijk koele 


| à , 
n . . gebied, onafhankelijk van de 

Ee hoz grootte van het gebied met de 

EN dichtheidsperturbatie op t= 0 
‚tegel , . e Î e se 5 Ei … a TP 
ngen ‚ s alleen afhankelijk ven DX en DZ, 
‚ Or OE ‚ Dit is een numeriek effekt, dat 
ì EE , de afmetingen van het uiteinde 
Ce sent rua; | | tijk koele gebied kunstmatig 
x=richting + J zonsoppervlak __ etabiliseert,. 


Falveert men DX en DZ, dan wordt het koele gebied corresponderend kleiner. 


d, e ziet men hoe de enelheid van de inatroming gericht 13. 
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Het pincheffekt wordt ok gedemonstreerd door de fig. 9, Hier is uitgegaan 
wan eeu zich over nog meer punten uitstrekkende dichtheidsperturbatie op 

ts 0. Ook hier blijven uiteindelijk 9 punten over waar da dichtheid hoog 
en de temperatuur laag ia t.o.v. de omgeving. Ia fig. Sa zien wa da dicht” 


a 


heideprofielea în de zerichting op verschillende cijdstippen. Hetzelfde Âe 


in fig. 9d gedaan voor de sneïheid. Fig. Sb geeft weer hoe da breedte van 


het afkoelende gebied afneemt als functie van da tijd. Hoe de temparatuur 


in de verschillende punten verloopt als functie van de tijd ie te zien Îa 
Lig. 3e. | Oe | | 
De en hiervan is als volgt: 

Beschouw de cellen, gevuld met koel gas, die grenzen aan het hete g42. 


mi 


Alle massa dis de condensatie binnenstroomt als gevolg van de afkoeling van 
da hele condensatie, stroomt eerst door deze cellen. Hier koelt het ges 
voor een groct deel af en vervolgens wordt het gas tot de temperatuur Ân 
het centrum van de condensatie efgekaeld in de meer naar binnen gelegen 
cellen van de condensatie. M.a.w, de randcellen worden heter door de 
centrale cellen, Hun stralingsverliezen zijn daardoor geringer (natuurlijk 
alleen voor temperaturen waarvoor de uitstraling afneemt bij toenemende 
vemperatuur, zie fig. Î) en zij koslen dus langzamer af dan de meer can 
trale cellen, Povendien daalt de dichtheid in de randcellen plotseling 
sterk als het proces niet meer isobarisch verloopt gedurende de zaalis 


fase van de afkoeling. Icmers, de randcellen zijn de eerste die in aap” 


merking komen om hum nog wel snel koelende buurcelien (die meer naar heft 
veutvum van de condensatie liggen) van massa te voorzien. Ds vandceelien 


worden vu snel verhit omdat: 


a} hun dichtheid lager is geworden en daar ook hun warmtecapaciteit is 
afgenomen; | | 

bì er veel coupvessiewarmte vrijkogt en het naar binnen stromende as 
heet Íe3 

2} zowel de afname van de dichtheid als de toename van de temperatieid de 
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uitstraling verminderen. 
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Warm tegeleiding 
EAT. End ie hen de denn delende 


Zoals we hebben gezien wordt het effekt van dervan leide wanndiltkhant 


met dat ven de afkoeling t.g.v. de storing in de dichtheid, als 
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Lost men deze vergelijking op voor Kk, (zie blz. à) dan windt mea 


f {Te 3 Ja 
K, as ai 


| ” 3 
ofwel met LL 1,5 DE 3509 kut , O, = 0,25 * 10° deeltjes/cm’ 3 


108 ° 
Te iC” Ks 
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Dit zou betekenen dat de warntegeleidingscoëfficiënt ongeveer 1069 kear za 
klein gekozen moet worden als de werkelijke wermtegeleidingssoëffiniënt ven 
de corona (bij afwezigheid van een magneetveld bij Te 10° K) ow voor het 
gebied met de gekozen afmeting {3 DZ X 2 DA) ex voor te zorgen det het 
thermisch instabiel wordt, Dit is echter niet helemaal juist, zoals blijkt 
aropvolgende numerieke resultaten. 
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Laax deze temmeratuurput als gevolg van warmtegeleidiag cvergaan in San 
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minder diepe en bredere temperatuurput, bijv. dubbel zo breed (zie fig. Â Ò 
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Yoe verandert de toraie uitstraling bij overgang van situatie TJ naar situ= 


atie IÌ (onder inviced van warmtegeleiding) 7 
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a 80% ela Om 
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e= 1,000 als a = 


bis Te i0f K geldt ongeveer ae =l : als den iets kleiner ís dan 

dan is a ®=3 , Het resultaat is echter niet aterk afhankelijk 
van de waarde van a . We zien dat de warmtegeleiding een herverdeling ven 
de wasgs cet gevolg beeft, die een wat kleiner netto energieverlies tct ge” 


nie d he ce Be 
vais heeft, 


dln als gevolg wan warmtegeleiding met een factor 2 verminderd. Cp den 
duur, als de temperatuurgradiènt voldoende verminderd is, kan de afkoeling 


door stralingsverlies het uiteindelijk toch van de warmtegeleiding Winnen. 


Gj 


Dit. alles wordt gestaafd door enkele numerieke resultaten. 

In Zig. la ís uitgezet de temperatuur van het centrale punt als functie 
wan de tiid voor de gevallen 1) zonder warwtegcleiding en 2) met warmtege” 
leiding; wacntegeleidingscoëfficiënt 0,0} van de warmtegeleidingscoëffici= 


eaf. van de corona zonder meagneetveid (bij … 10° K) . 


le zien det in het tweede geval de temperatuur eerst langzaam een weinig 
oploopt, dan vrijwel constant biijft en tenslotte snel afneemt, Als een= 
maal de theraische instabiliteit goed inzet, verloopt alles geheel analoog 
aan het geval zonder warmtegeleiding. Het voornaamste gevolg van de warmte” 
geleiding is hier den ook dat het veel langer duurt alvorens de thermische 
instabiliteit inzet. In fig. Oc, d, e zijn op enkele gehen 
fielen, teuperatuurprofielen en netto energieverlies (per cm° } profielen 
weergegeven, benevens het verloop van de druk en de Menebeid van het cen 


trale puut als functie van de tijd (fig. 1Ob}. 
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Daar op t = Ô de temperatuurgradient in de xerichting een faktor 
In Ei ® Le € . 2 e & e Ei ie 7 
B” 35 kleiner is dan in de zerichting, vindt in deze richting nauwa= 


lijks een vervlakking van het GRRUSEAENER OTR plaats, 


Laat men hetzelfde proces verlopen bij een ongeveer JX zo hoge dichtheid, 


waar met dezelfde warntegeleidingsccëfficiënt en dezeïfds afmetingen van 


bet koelende gebied, dan merkt men bij 10% storing wan de dichtheid ín 


Dd 
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het geheel niets van de warmtegeleiding: de thermische instabiliteit zat 


dir ect in. 


Verschillende stralingsfuncties 


in fig. 31 zien we de temperatuur in het centrale puat uitgezet als func 
cie van de tijd, echter nu uitgerekend met de Mo Whirter stralingsfunctie 
voor verder hetzelf de probleem ala het voorgsande zender warmtegeleiding. 
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De Me Wairter stralingsfunctie, die uitgaat van zoane=sbondanties (zie 
bijw. 10) wijkt nogal af van de Cox & Tucker stralingsfunctie (zie fig. 


1). We noemen bier de voornaamste rh nd 


“8 5 « ml cW -e e | Pi C & T 
1 Bij T= 10° is de waarde van £° ongeveer de helft van £ : 


Sovendien neemt de Me Whirter funcecie veel jan ngzamer toe bij afname van 
de temperatuur. Dit heeft tot gevolg een veel trager inzetten van de 
instabiliteit, 


2) Het waximum van de Me Whirter functie ligr. veal Jager dan dat van de 


c3 


ox & Tucker Functie, Bovendien is het maximum van de Me Wiieter func= 


tie verschoven naar een lagere temperatuur. 


Go 


) Beide stralingsfuncties vertonen een piek hij ongeveer 20 OCO K aie 
gevolg ven Ho , Deze piekwaarde is woor de C&T functie echter veel 1a- 
ger dan voor de Me Whirter functie en bovendien valt da piek van ds 


Me Whirter functie bij een iets lagere temperatuur daa wan de CAT funce 


EL 


Zoaïs reeds eerder vermeld wordt de plek van de uitstraling per volume 
eenheid in bet geval van een isobarisch proces bereikt bij het doorlopen 
van deze plek rond 20 000 K van de stralingsfunctie, Dit is weergegeven 
in fig. i2a, In fig. 12b is hetzelfde weergegeven voor de uitstraling van 
bet centrale punt, zoals berekend met SHASTA. We zien dat het verloop gla= 
baal hetzelfde îs als dat van het zuiver isobarische geval, alleen Be 
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Bij temperaturen lager dan 26 000 ‚jn de resultaten niet meer te ver 
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trouwen omdat er dan een aanzienlijke hoeveciheid neutraal waterstof ge= 
vormd gaat worden { N(R) / ne > 1,09 woor T e= 2. 10° K} en omdat 
de dichtheid bij die teaperatuur ongeveer met een factor 50 is opgelopen. 
Dit veroorzaakt een verstoring van het proces als gevolg van recombinatie 
an het optisch dik worden. Bovendien woet (zie later) de materie van de 
condensatie deor een magneetveld worden ondersteund, hetgeen echter voor 


neutrale deeltjes niet weer opgaat, 


Masgef lux 
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In hef tweedirensionalie geval wordt er massa aangevoerá vanuit twee riche 


tingen, de x= en de z=richting (zie fig. op blz, 14 en de figuur met vec 


torpijlties), Zij M de massaflux vanuit de x=richting en M, vanuit 


de zerichting, gedefinieerd alg: 
3 
M e= 2x Ì {Pp Vv). & DZ 
x : KR 3 
iel a 
3 
M= 2x E (Pp vo), * DX 
ij ie} ì 
zedat 


4 LCP vd, 


MDK ' E (pv), * 
k Hd Kl 
M 
Het verloop van rin is te zien in fig. 13 in het geval van de Mc Whirter 
ee ee Dn 
stralingsfunetie en were 0,3 . 
de EN 


% Ld Rl had 
Ëert na € = } 18 


Li 


= 0,t , loopt vervolgens op tot ongeveer 6,25 
Le0 Ís de devkgradient in de smalle (= z-) richting groter dan in de 
brede (= x=) richting, immers hetzelfde drukverschil Ap vindt plaats ie 
de gerichting over DE en in de x=richting over DX , met DZ < DX , zodat 


aanvankelijk wv, groter wordt dan vw, . 





3 Ee BE 
Echter hierna gaan de traagheidstermen Vv 57 en Ov, ae ri in de bewer 
Ee 
gingsvergelijkingen een vol spelen: daar in het centrale punt zowel zn 


ik 
als v_ pul moeten zijn is er een gradient van Vn in de z=richting en 
an 


od 


DW a Ea fs A er . Ro « 6 DE ami e 
baj 500 CO9 £ en valt dan weer af bij nog lagere temperaturen, Direst mn: 


mtr 


van v, in de zerichting (zie snelheidsprofielen van Î dim, berekeningen). 


Geze term is voor de x-richting veel kleiner dan voor de zerichting, zodat 
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drukgradienten in beide Fichtingen zijn ongeveer hetzelfde geworden. Dit 


| PN uri À 30% PT ER Get : er bh nf 
„äatste wordt anders als het proces niet meer isobarisch verloopt: dan 


A 
o 
À SL Le’ er DE - é KE, dn 
Be cer dan ze en wv, dus ook M, beginnen toe te nemen 
Lee. OV EESD. M_. 


kee oe De #5. wim dn ae ,  « ie , 6 : 3 En ne ps 5 
aar het proces gelijk veel aneller verloopt bij gebruik van de CAT atra= 
: 4 nde en Ee Es k D $ e 9 
ingsfunctie is u” daar vrijwel constant, ongeveer 0,35 
. ab A " 


Rn, | 
kml 


& 


EER OR Oe RSE 


EET A A ad in nnen hak RETE RT Vk TS nn lk id  id 





Een an alvtisch: modal voor een thermische instabiliteit ( die plaatsvindt in 
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We zullen hier twee wijzen oaderscheiden, Ten eerste cen iscdensiech 
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iteit en ten tweede aen isobarisch verlopen” 
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Op tes Q wordt in een volume V sen dichtheidsstoriug van 102 aange= 
bracht. We nemen aan dat het proces zo snel verloopt, dat de dichtheid 
in V en erbuiten niet noemenswaardig verandert aîs gevolg van de af- 
koeling. Hoewel we op ts= 9 uitgaan van drukevenwicht, is er gedu= 
zende her proces geen sprake meer van drukevenwicht, Nadat de afkoeling 
is gestept (d.w.z. na een tijd Te ) zal het drukeveawicht zich wervolgens 


d rt, A . {L is karakteristieke afmeting V en 
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0, = dichtheid in V (onafhankelijk van de tijd). 


Dm dichtheid buiten WV (onafhankelijk van de cijd), 
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en RAD (p,T) = netto energieverlies per sec Dar ca 
ais gevolg van uitstraling en wecha- 


nische verhitting 
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Dan laat de energiebehoudvergelijking zich schrijven als 
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De beschrijving van het isobarische geval Îs zeer veel lastiger. 

De ee aannamg îs julst ale NE ‚. Het afkoslen vaa V 
heaft nu een grote maasaflux tot gevolg, die de druk ondanks de 
temperatuurdaïing constant haudt. Er wordt au dus hete materie wan 
buiten V aangevoerd, Deze materie moet worden afgekoeld en gecom” 
primeerd. De energie=afvoer vindt alleen plaats iu de vorm van uik 
straling. 

We maken de volgende veronderstellinges: 


%, KT Se U 
„ de Sruk p= En blijft zowel binnen als buiten constant gedu 


rende de vorming van de condensatie. 
… de dichtheid o{t} en de temperatuur Tt} zijn uaiform over Y, 


Buiten V. zijn de dichtheid P, en de temperatuur Te onafhenke” 
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lijk ven de tijd, 


We beschouwen nu een bepaalde hoeveelheid gas buiten V , die een von 
lume VV  , inneemt, Als dit gas V binnenstroomt worát het plate 
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seling sterk gecomprimeerd en vervolgens afgekoeld, De compressle 


heeft tet gevolg dat er een anergie Aw = pAV vrijkomt, waarin AV 


de volumeverandering is van het gas dat we beschouwen als gevolg van 
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ionen) in V inneemt. Daarnaast is de thermische aneseiedahoaa VAL 


de beschouwde hoeveelheid gas veranderd en wel met een bedrag 
AE = Er 4 B ee md AV, 
Vor had, Pasie vel 6 ae 
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Volgens de eerste Hoefdwet van de the ermodynanica gzeidt nu: 


waarin ÁQ de toegevoerde warmte is. 
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Als we voor AV kiezen de verandering van het volume die het gas dat 


in Ì sec V binneastroont ondergaat, dan is: 
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Daamsast moet het gas dat zich in V bevond voordat het volume AV 


erbij kwam worden gecomprimeerd en afgekoeld, Voor dit gas bedraagt 
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Dit moet worden afgevoerd in de vorm van straling (met daarvan afge- 


trokken de mechanieche verhitting). 
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net weer BAD (p,T) = E{T) p° — £CT) ai p° 
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Nemen we de serete orde vartatie van g(x) mee, d.w.z. 
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Voor een 102 storing op Teu 0 vinden we hieruit (tT << 1) 
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Hoe verloopt de {isobarische) afkoeling van een element dat men volgt? 
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u. Hoe beïnvloedt warmtegeleiding het verloop van de afkoeling in geval 
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De evergievergelijking wordt: 


V . RAD (o,T} = kv. Sen TS. 0 ; 


aarin CQ = opperviak van V . 
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We beschouwen mogmaals in het kort een aantal aspecten die van belang kuna” 


sen zijn woor het verloop van de ontwikkeling van een rustende protuberans. 


Ì. 
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6. 


=l 


Opdat de thermische instabiliteit in kan zetten, is een lokale dicht= 
heidsverhoging noodzakelijk. Deze zal in het algemeen isobarisch zijn 
vanwege de hoge geluidssnelheid en/of de hoge Alfvénsnelheid, | 

In een neutrale laag kan da dichtheid tot veel hogere waarden oplopen 
dan daarbuiten als gevolg van drukevenwicht en dit kan resulteren in een 


veel hogere uitstraling aldaar, 


De sterke warmtegeleiding in de corona zal de afkoeling verhinderen, 


er 


enzij de warmtestroom wordt geblokkeerd, of het afkoelende gebied zeer 
groot is, (Zie biz, 5, een gebied met een grote dichtheidsperturbatie 
over een oppervlak van ongeveer 101? km? is vereist). Dit laatste ia 
enwaarschijnlijk. Oek hier biedt de ideale neutrale laag een oplossing, 
daar het magneetveld in staat Îs de warmtegeleidingsco&fficient in de 
richting loodrecht op de veldiijnen verwaarloosbaar klein te maken, Pro= 
blemen leveren natuurlijk de onder- en bovenkant van de neutrale laag, 
daar hier wel warmte naar binnen kan stromen. Het ie echter niet ondenk= 
baar dat er meer ingewikkelde magneetvelden kunnen worden geconstruserd, 


die een klein gebied voliedig kunnen afschermen, 


Een mogeiijkheid hiertoe wordt bijvoorbeeld geboden door de tearing-mode 
instabiliteit (zie Kuperus, M, and Tendberg-Hanssen, E‚, 1967 en Appen= 
dix, biz. A3), die aanieiding kan geven tat een ge esicten veldconfigura= 
tie, 


Een eindig elektrisch geleidingsvermogen kan bovendien een sterke anni- 


hilatie van magneetvelden tut gevolg hebben in een klein gebied rond de 


neutrale laag, hetgeen aldaar tot aanzieniijke vertraging van de afkos= 


Ting zou kunnen leiden. 


De ideale (verticale) neutrale laag kan de zich verdichtende materie 
niet dragen, mits ket magneetveld ook in staat is de protuberans te on= 


derstennen (zie Appendix, biz. A2 en 9). 


De veetpuanten van de veldiiinen zitten vast gan de fotosfeer (line= 
tying), hetgeen de vorming van condensaties near ouderen toe bemoei 


Likt. 
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Naar beven toe neemt de massadichtheid af, hetgeen eveneens de vorming 


van condenssties zal verhinderen op grote hoogte. 





28, 


8, De mogelijkheid tot condensatievorming hangt sterk af van dae abondaaties 
van enkele in de corona schaarse elementen, van de temperatuur waarbij 
de thermische instabiliteit moet gaan inzakken en van de (lokale) mecha 
nische verhitting, 


Vooral de waarde en de gradient van de st-ralingstunctie bij de begiaten= 
peratuur zijn van groot belang voor de bepaling van het al dan niet in 
etten van de thermische instabiliteit, Voor de Cox & Tucker kromme zien 

we (z ie P.P, Cox aad W.H. Tucker, 1969) dat dit (bij io% &) vooral worde 


bepaald door de abondartie wan zuurstof en neon. 


9, De eindfase van het afkoelproces wordt baliemerkt door recombinatie en 


bet optisch dik worden van de protuberansmaterie. 


De meeste complicaties die zich zouden kunnen voordoen bij de vormieg van 
aen rustende protvberans zijn we hier uit de weg gegaan. Niettemin : zin 
enkele van de bovengenoemde punten aan de orde geweest. We hebben gezien 
âzt een kieinere warwtegeleidingscoëfficiënt als die van Spitzer voor de 
koroná de afkoeling niet hoeft te verhinderen, doch wel aanzienlijk kan 
vertragen, Verder bleek uit de berekeningen dat de condensatie de neiging 


e 


vertoont om voortdurend smaîler te worden ala gevolg van een achterblijven 


net koeten {vergeleken met hei centrum) van de buitenste gebieden van de 


had 


condensaties, Verder zagen we ook de invloed van de abondanties van de 


ï 


ekementen op het koeïprcces, zoals geillustreerd aan de hand van de Cox & 


fueker= en de Me Whirter-funotie, 


De berekeningen die hier zijn gemaakt en net name de tweedimensionale be 
rekeningen, zijn uitgevoerd met een zeer gering scheidend vermogen. Een 
beter scheidend vermegen vereist meer punten en een kleinere tijdstap 
(stabiliteitsconditie!), hetgeen tot buitensporig lange rekentijden aan” 
beiding geeft (zie Appendix), Tech ís een goed scheidend ve rmogen van 
groot beiang voor bijvoorbeeld de beshesiodns van een ingewikkeld mag= 
neefveld, diffusieprocessen, de tearing=mode instabiliceit: waterie in een 


zwaartekrachtveld ate, 
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a} De continuiteitsvergelijking 
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e} lgestandsvergelijk: 
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5 = magueetvald 
g = versnelling van de zwaartekracht 
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p nn: 54) 
2m inwendige energie per massa=eenheid = Ee (y= 5/3, 
Kk = wermtegeleidings svermogen ideaal gas} 
8 kod Ed ud €. 
ns resistiviteit = Tra Pet JS = electrisch geleidingevarmogen 
Î } 
e= lichtsnelheid 

H = gemiddelde massa per deeltje 

MECH . 3 » NE: s Dan 4 men 
Ö = mechanische verhitting per volumereenheid per tijds= 

… eenheid 

AAD 
d, = uitgestraalde energie per veolumeveenheid per tijde 


eenheid, 
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De viekositeit van het gas ia bij deze vergelijkingen verwaarloosd, 
Bij de ae berekeningen ís de invloed van de zwaartekracht niet mae= 


gevomen (dus 2 zen Ö), evenals de magnetische diffusie (n = 0}. 
De zwaartekracht: 


Om na te gaan of de invloed van de zwaartekracht van belang is, beschouwen 
we een geschematiseerde condensatie: cen vertica ale cylinder, straal R, 


iengte L (zie fig.) 





dan 


Zij wv. de snelheid waarmee materie zijdelings wordt aangevoerd naar de 


condensatie toe, en rn de snelheid waarmee de materie onder invloed 
van de zwaartekracht naar beneden stroomt. b. is de dichtheid in de ey= 
binder, Pp, de dichtheid buiten de cylinder. Opdat er meer of evenveel 
materie zijdelings naar binnen stroout, als er door de zwaartekracht cat= 


rrokken wordt, moet gelden: 





re Es : 
Zn R.L. vw. > Ov „ERS f9 — 0 
in ° Pe gut He 09, 
Ëwei 
A 7, 
$ x Vv 
He out 


zal meximaal gelijk zijn aan de geluidssnelheid, terwijl Von Snes 
veer gelijk zal zijn asu de vrije valsnelheid van materie die over een af- 


De 
stand # L versneld is,onder invloed van de zwaartekracht g 


max PE 
Dus Vv “y_ en v 4 / LG 
in s Gut v ° 
Vullen we dit in voor een concreet geval, bijv, E= zE b. 
e en ’ sE AA 
vO 1,5 % 167 cm/sec s L= 10°? em Ss 2 = 2,5 # 19% em/esec? « dan 
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O.wez. als de dichtheid in de condensatie groter wordt dan 20 keer de 
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je aan de wensten geen stabiele situatie kunnen ontstaan: de 


condensatie lekt leeg onder invloed van de zwaartekracht. Ala bijv. 
) 
We, & v wordt deze conditie: aj Î.2 , zodat dan zelfs een 
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dichtheidetoevame van meer dan 204 niet meer ín stand kan worden gehouden, 
De zwaartekracht speelt dus een belangrijke rol, zodat een goede theorie 
over protuberansen een mechanisme ter ondersteuning van de peers zal 


moeten bevatten (zie o.a, DP), 


Oueindig geleidingsvaermogen 
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Magnetische diffusie wordt beachre en door 
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Da karakteristieke tijd voor diffusie is dus 


22 
Z * e pn . a 
En u men Met L, « karakteristieke afmeting waarover B aanzien 


lijk verandert = 3500 km. 


mt & Beed 3 
Ha is nn 6,5 w tai? 


Dec met HA 19,3 , dus bij Te 10e xs 


ne 5 
n & 1,3% 10" e.m.u. 
Hieruit volgt voor magnetische diffusie de karakteristieke tijd 


cm 101% sec 


fn 


Vergeleken met de tijdschaal waarop de vorming van een protuberans plaats” 
vindt, = 10 sec, is de maegnetieche diffusie dus niet van belang zeleng 


er gesu sprake is van een avomaal geleidingsvermogen. 


Een andere instabiliteit die kan’ optreden als beval van een n £93 ás 


de tearing-mode instabiliteit, die een aanzienlijk kortere groeitijd 
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heefi: 





| Es : 
Tr. 7 VTT, met T‚ = Alfvéntijd = edt) Ld Wai e 
f } Î b N 
LT & Â Aje) 
À Er pad hd ET epe 5 5 3 Ù UE 3 . en à ter ; f 7 / á ’À ä ke 
Âls B 5 Gauss , n, = 0,35 x 10° cm ° 5 v‚ “6% 10 cm/sec , dan 18 
TG sea EL Tw TE 2,4 X 10° sec. 
EN kelde 
Dus cok re M. is groot vergeïeken wet de karakteristieke afkceltijd. Pas 


P ani 


bij afrmetin van het koelende gebied van enkele honderden kilometers kan 





À van de erde van grootte van de'groeitijd van de thermische instabi= 
titeit worden. De beschrijving van dergelijke kleine gebieden vereist ech= 


n hoog scheidend vermogen (d.w.z, kleine staplengte) en, vanwegs da 


F% 
5 
a 
® 
® 


a 


stabiliteitsconditie van het numerieke schema, een kleine tijdstap. De be- 


Ld 


nodigde rekencijd is fa de praktijk niet te realiseren (nie later); zie 
„ak 8. Si ASTA. | 


PE 


Aije op te lozsen vergelijkingen hebben (in één dimensie) de gedaante 


3E 8 

zen en gen (Pay) e= g 

de dx ) x 
Hierin kaa PF geen van de volgende grootheden zijn: @ (eontinuiteitsver= 
gelijking). Pv of Pv, (bewegingsvergelijking in resp. x° Of gertech= 


x z 
Ling), 5 ot B (vergelijking voor x= of z=component van het magneetveld) 


led 


cf pe (de energievergelijking). 


*. 
ed id Son, 
a representeert de overige termen, die naast termen wan de vordt xp € 


3 
5x (Fv) in de vergelijking voorkomen, 


Een der mogelijke numecieke sctiema's waatmee deze vergelijkingen door iten 
veren kunnen worden opgelost is het Lax-Wendroff schema. Dit schepa vol= 
doet echter niet als er steile gradiënten in de te konvekteren grootheden 
voorkoren. Deze steile gradiënten worden namelijk snel vervlakt als zevolg 
van de sterke numerieke diffusie, die een eigenschap is van het Lax-Wend= 


Toff schema en ven vele andere numerieke schena's 


Tijdens het opzrede n van de thermische instabiliteit en daarna is er apra= 
ur, hetgeen niet verwon= 
derlijk is. Immers, het eindsrodukt is een koel akkak var hoge dichtheid, 
3 ze b d Ae es T$ h es , Kat „% a à de dt OP VEN BPN ie 1 „nn Ed de 
ingebed in een hete corona, waarbij de dikte van de overgangslaag van 


ie condensatie zaer veel kleiner is dan de dikte 


« « 


Het is dus gewenst een numeriek schema te gebruiken dat wel deze steile 


gradiënten aankan. Ren mogelijkheid is het schema SHASTA (Strong and 
Smooth Transport Algori.thm), dat gebruik maakt van een techniek die bekend 


steat als FCT (Flux Corrected ideen het kort komt het er op zeer 
dat eea tijdstip uu uit de volgende gedeelten bestaat: 
1} konvekteer E bijv. met het Lax-Wendroff schema, 

Tijdens deze stap treedt sterke diffucie op, die ongewenst is, 

Daarom wordt er een tweede stap toegevoegd? 
2} Anti-diffusiestap. | 


Pijdens deze stap wordt de diffusie uit Ì) ongelaar gemaakt, echter zo= 





ú 
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DS | 


Ba 


danig dar er geen extrema gegenereerd of versterkt worden, hetgeen 


wordt bereikt door corzesties uit te voeren op de massaflux (FCT). 


Tweede-orde nauwkeurigheid wordt bereikt door de brontermen uit te rekenen 
iidetip t + Ì DT (d.w.z, voer da konvektietijdstap uit met als tijd= 
stan j UT en de brontermen berekend op tijdstip t) en deze te gebruiken 


had 


je 


hij de kouvaktie met tijdstap DT , Voor meer details zie 2, 3 en 


Da ViTennthede heeft als nadeel dat her ook een nieuw miaimum of Max LRAD 


Zat freiech veel íe, sfkapt als heat zich over woinder dan 2 punten uit= 


‚ hetgeen vooral vervelend is als de dichte, koele con= 
denaatie, za de thermische pinch, zich nog slechts over enkele punten uit 
strekt, Wet ovlossen van de vergelijkingen in twee dimensies verloopt ge- 


hee} acaloog door toevoeging van een vergelijking voor de andere dane 


voe vn fen hebben we aangenomen dat alle termen van audere gedaante als ze 


A 
en Ey Aaen a 2 Là a 
ok a (Fy) wearwerkt zijn in de bronterm 8, © 
si & eT 
En Lt, TEA 5 s B o @ 
De Aiffusietern see ÍK se) woor cel ij kan mee met een expliciet schema 


9% 0 dx 





{ wr ef X Ez ae T 
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; : : oep? iel 1 Ee EDITED ERIK IDEEEN 
p { es úÂ Ì A+ ä be Éx Pi $ \ DX 
eg & — ‘ ” 
Ke) C 3 DE 


met 5 4 (K.K, en DX = HK, 7 Ke 
hadt Er k $ Ü Ei “Le? + hei & EN i 2 


Indiea de invloed van de brontermen zeer klein is vergeleken met het ef- 


Fekt vaa de konvektie, is de stabiliteitsconditie: 


eg 
hed 
jeej 
EN 
et 
tas 
te} 


ee se Ominimwim ven (DT DT ler Ò 
MAX gt et EM 


| | Dx DZ 
DT Bà Gi m ur nd e DT SER ë Kr nd EET KD 


Vs) ax i 


Lr {DX D4 jd . 
DI, mT a» nin DEDEL met Vp u ° Fast-mode snelheid = 
“bh 4% F.M. Kd 2 ananda ttmnndbnnnnde 

aal ve + va 


geluidesne elheid = DR. 


EE 


KA 





„ 8 2 5 
Alfvénsnelheid = —mr 
| v &mp 


< 
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Courant getal. 


Le, 
# 


ie alens, ESE ak Tk EP ir va en ET Ted dc nd 


Te rens GEA ke wee Sonde ene sene en 


SE dat 


5 


w zeker te zijn van stabiliteit kan men het best & S 0,5 kiezen. Het 


optimum, wat betreft nauwkeurigheid en rekentijd, ligt bij ongeveer 


ú 


a CZ , Als men de bronterm voor de warmtediffusie van blz, A5 in da 


vergelijkingen opneemt, komt daar nog bij de stabiliteitsconditie: 


pr < 0,5 CUL 


KMAX 


à 


Vooral bij. hoog scheidend vermogen (DX kiein) en een grote warmtegalei= 


dingscoêfficiënt (,,, groot) kan dit zeer grote, ongewenste beperkingen 
oplegsen aan de tijdstap. 


De stabiliteit wordt ook voor minder sterk beperkte tijdstappen gewaar= 


borgd, als men de diffusieterm in een aparte diffusievergelijking atopt: 


Tr 3 {ar} T 
re Led se | AT $ „je s 
Ö 35 À Öx J & 


die via een impliciet schema (inversie van een tridiagonaalmatrix) wordt 


3 TE 5 n 
gogeiost, De is een bronterm die eventuele warmtebronnen (cf -putten) 


re 


ke : À > @ © 4e n @ 
in cekening brengt. De totale bronterm s wordt daaxbij gelijkelijk 


Sea Ean a NE: 
over Er e€a zerichting verdeeld : 8 OT mies 


T 14 Sei 
De bronterm s kan ook In de brontermen van de konvektievergelijking 
worden gestopt, echter is het in verband met dae stabiliteit beter deze 


hrontermen in de diffusievergelijking te stoppen, 


T „2 0 9 nj s Kd 
Äts de bronterw a zeer groot wordt (bijv, in het geval van de thermi= 


ge 


sche lustabiliteit als het netto=energieverlies per volumeveenheid ín ‘een 


ee 


tijdstap van dezelfde orde van grootce wordt als de thermische energie= 


) 
inbeud per volwaeeenheid) wordt de stabiliteitsconditie 
DT << ram mek ĳ << 1, 


(In het als voorbeeld gegeven geval is at = (QT =@Q 7) en 
er 
pa == =p} 
2 
De konvektiestap wordt uitgevoerd na de diffusiestap. Dit leidt tijdens 


de konvektie niet tot instabiliceit zolang DT zo is gekozen dat het 


effekt van de diffusie kle 


Ee 


Da gebruikte computer ís eem CDC 6500 (ACCU, Utrecht). 
in onderstaande tabel is de CPU-time per tijdstap voor deze computer voor 


enkele roosters aangegeven: 


veel kleiner worden dan de DT ax Var blz. A5, hetgeen natuurlijk wel af 


hankelijk is van de gekozen DX , a en de optredende wv 


cuschrijving rooster __ CHU-cime per tijdstip 
OPLIEPEN BETTEN BED EROETRIKT RCREEE EDEL RTD ON RATTEN CRIME 
iD: 6) punten | 0.3 sec 
DD: 100 punten 0.5 sec 
2D: 2Î Xx 15 punten 2.2 sec 
2D: 63 X 43 punten 17.0 sec 


De totale bensdigde rekentijd hangt nu af van: 


1) de ataplangte DX (of sok DZ voor 2D-problemen). 


« 


DX bepaalt niet alleen het aantal benodigde punten, doch ook da SEOOREG 


van de tijdstap, zoals bepaald door de stabiliceitsc onditie 


2) het verloop van de uitstraling als functie van de tijd: 
Als het netto=energieverlies als gevolg van sterke uitstraling namelijk 
erg groet wordt, moet DT worden aangepast (zie onderaan vorige bladzij= 


de). De tijdstap wordt dan DT 4D . 


DE an ken gedurende de sneile afkoelfase (10 009 < T < 60 600} zeer 


Apen 
MAX © 
je tataal benodigde rekentijd voor het doorrekenen van het condensatieprer 
s hangt dus sterk af van de begincondities, met name de begindichtheid, 
de begintemperatuur, de storing in de dichtheid op t = © en bec verloor 


san de stralingsfuactie van de temperatuur, 


Voorbeelden rekentijd: 


a) Eëndimegsiornale berekening (zie figuren 5 en & voor gegevens) met 
x= 0.0 5 a = 0,5 … 
Er waren ongeveer 14 000 tijdstappen nodig, d.w.z. ongeveer ÌÀ uur 


vekeutijd,. Hierbij is het scheide nd vermogen groot. 


bh) Tweedimensionaal 2Ì X 15 rooster (gegevens, zie fig. 10, zonder 


9) 
warmtegeleiding, Cox à Tucker stralingsfunctie}; rv * 0,05 3; am 0,3 . 


1 
Afkoeling van 10° tot 10* K vereiste ongeveer 3500 cijdstappen 


gn 
Nae: 
Dd 
te 
jet 
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au met Mc Whirter stralingsfunctie, 


+ 8500 tijdstappen —… 5 uur 


Opmerking: Zou men het probleem van a) met een 63 X 43 wooster wijlen be 
schrijven, dan moet de tijdstap ongeveer drie keer ze klein war” 


t 


den gekozen om sta 


& 


biliteit te houden (immers: DZ is ongeveer 
drie keer zo klein) en er zijn ougeveer tien keer zoveel punten, 

„ Ee © « ft „IX 
d.w.Z. er is ongeveer dertig keer ze grote rekentijd (+ 60 VUL e j 


nodig. 





Gegevens be behorende bij de. de figuren: 


dee 


’ 
| 
Ì 
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Figuur Îz 


Deze figuur toont de gebruikte stralingsfuncties (T). Zij zijn beide 


genomeerd op Ì voor een temperatuur van 10° K en tonen de temperatuur” 


__… afhankelijkheid van de uitstraling van eea heet optisch dun plasma voor” 


a 
RRC: 


ee « 


resp. kosmische abondanties (Cox en Tucker} en zonneabondaaties (McWh: 


sek 


De uitstraling zoals gegeven door de Cox & Tucker=functie ligt hijme over 
het gehele temperatuurbereik boven die gegeven door de MrWhnirter-funsctie. 
De werkelijke waarde van de Cox & Tucker functie hij 10° XK is bovendien 


kf  « 


. | bijna dubbel zo groot als die van de McWhirter=funetie baj io E, 


el 


à& + 
uur 4: 


armnamewern Â PART 2 BAP RDS" EED 


| | Thermische instabiliteit bij hoge dichtheid (d.w.z, bij een dichtheid die 


ker ie dan die van de corone) Cox & Tucker stralingsfunctie. 


ws 2,5 km; eenheid van RAD e= 339 erg/cm’ sec). 


X 
Positie O = centrums positie 40 = rand rekengehied (rooster); symmetrie 


honities G t/m 14: 5% hogere dichtheid dan de punten daar” 
buiten 
| 10° 
: temperatuur = T, = aen Áo 
4 . k, iS C&T 
mechanische verhitting = £° OA 
posities 15 t/m 4Osdichtheid = p = 10 ii g/e 
| ne Bo: ie 
| temperatuur = T_ = 10° £ 
fl dn 
| ON dd C&T 2 
| mechanische verhitting = £ CE) JE 
| Ao tewperatuur en netto-energieverlies per volumeeen aid ala fupctie 
| van de tijd. | 
| b. terperstuurprotielen, 
C. snelheidsprofielen (een positieve snelheid betekent: beweging naar 


d. ennen Jen. 


1 Le 
e. drukprofielen (er is geen drukevenwicht vanwege de gekozen hoge 


ER ee A RL lk CH in a dk hik 2 it ik a an: he oi 


TREEDT rn B en 


2 8 z 


dichtheid; de karakteristieke afhankeiijkheid is veel kleiner dan 


de tijd waarin een geluidsgolf vanuit het centrum naar de rand reist). 


8 reds 
& 


netto energieverlies (per sec per cu’) profielen (RAD = uitstraling — 
mechanische verhitting); op t = 0 is RAD = } in het gebied met de 


dichchei idsperturbacies. 


med van een temperstuurput (950 000 X tegen omgeving 106 K) ale 
lin van renate Lag. 

tusevenwicht op & = 0. 
Geen atralingsatkoeling. 


ez * ap ê En ES Pa” ‚ Nun Ì 5 : Kl 
Dichtheid = 3,67 « A9 g/cu 


« bd 
a, erpsratuurprofieien op enkele tijdstippen. 
Gj; drukprefislen 


L, Svelheidprofielen (anelheid positiefs beweging van het centrum vandaan). 


Figuur dis 


MARS WP: Tie. AG DE ENT ARG ren S 


Fherwische instabiliteit + warmtegeleiding bij hoge dichtheid, 


Se, ,‚ ie, ze » 8 « % sd ®% K 3 é 
Cox & Tuekear strelingsfunctie. 


ka a va f ao vp sen 
be vand var 


keel 
“ 


het rekengebied en het centrum koelen kort va elkaar af doordat 


fe}: 


de warmtegeleiding het koelproces cok sterk in net gebied waar op tm 


geen 54 storing in de dichtheid was. 

, beed 3 Ë & ad er 2 ian ij - €: n es d 5 

Kr Ko 8 Kn FT 1Ô  ergeBes Ke warmtegeleldin a Lent. 
et EN 


Berekend m.b.v. een expliciet schema, gebruikt io de bronterm van SHASIA 


5 eN 


Hie appendix). 


ì 


Brukevenwicht op t = 0 (p = constant = 831,7 dvnefkm°). 


R&D ie per def. 1 cp t = 0 in het gebied met de 52 storing in de dichtheid. 
A. Temperatuur van het centrum (punt 0) en van de rand (punt 40) als 
functie van de tijd. Ook RAD voor punt 49 is geplot. 


b. tenperatuurprofielen op enkele tijdstippen. 
1 


en (snelheid positief: beweging naar het centrum), 


3 


d., nettovenergieverlies {pek cu° per sec)=profielea. 


ie 


‚ __drukprofielen. 
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Thermische instabiliteit zonder warmtegeleiding pb: 
dan in &a voorgaan nde gevallen. 

Cox & Tucker straliagsfunctie, 

Geen warmtegeleiding. 


DE = 10 km. 


drs 


vanten 51 t/m 53: 10 % verhoging van de dichtheid (f.Ì Pe 


puuten 54 t/m 55: lineaire afval van de dichtheid tot Pe 
punten 56 t/m 1015 dichtheid p, = 8,35 # 10 Es rd 
- De wechanische verhitting is overal 0.98 & EE, Joe zod 


cok de rand vog een weinig afkoelt 


Lagere dichthei id 


) 


a, temperatuur, dichtheid en Genen nea als functie van de tijd. 
END, , = temperatuur tu punt 5! 
TERE 3 a kemperatuur ju punt 53 
RAD, = netto snergieverlies vaa punt 51 {per sec per ca°) met als 
genheid RAD, (t = 0) =£ (T/l.1).(l.l p,)? = 0.98 ECT )p,P = 
… 2Î e e ee 
i.3 #30 3 erg/ (em sec}, 
KHO, “ dichtheid in puut 51 met ais eenheid RHO, (t = Oe 
bk, druk en aneiheid als fuuctie van de tijd. 
Pos 28 druk in puat Sis ep t = 0 ia Pe 51 overal 1. 
v(t} « maximale snelheid, d.w.z. de snelheid op tijdstip t waar de 


materie het snelst naar bet centrum beweegt. 
Ec, de verhouding vAv/ RAD ais functie van de tijd. 


mens 


Figuur 6: 


bamsdhtasdhans Heden nnn 


Kd 


Zie gegevens behorende bij figuur 5. 


. 


Ld 


Teuperatuur=, snelheid= en dichtheidprofielen op eakele tijdstippen. 


Figuur, EA: 


nek En 


A 


Thermische instabiliteit donder warmtegeleiding in een magnetisch 





Ed 


neutrale laag. 

Cox & Tucker stralingsfunctie. 
Punz 8 = „centrua, 

OÂ = 5 km. 


Eeginc cadities: drukevenwicht 


< 


hi an 
wasneetveld ; Bíj) BAK # ATAN (0.072 % r: 8D) 


1570796 
BMAX = 5 gauss 


EE: : d 5 fd 
Dichtheid ver van de neutralte laag: Ph 1.67 2 1Q=ek g/cx 
Dichtheid in de neutrale laag: Pe » 8,10 On 


RAD (t = Ö, neutrale laag) » 4.9 « 10 } evg/ (eu gend. 


RAD {t = Q, rand rekengebied) = 0 


we 


s ne & 8 
Begiatemperatuur = overal 16° KE. 


& 


De ploties tenen dichtheid, snelheid, temperatuur en magveetveld als furctie 
an de afstand tor de neutraje laag. 
De laatste plotjes van deze serie tonen profiel: an op t = HI,Î/ see 


berekend met een staplengte DX = 5 km en op t = 513.02, 
En 


ot t = 89 sec eu vervolgens DX = 2,5 


Di 
Merk op de smallere dichtheidspiek en de smallere tempers 


: % 
twaede geval, 
ki 2 en 


Tweedimensicnaal rooster, 15 x 15 punten. 
ge k Tucker stradingsfunchie, 
cen warmtegeleiding. 


Kr Rt % pees 
ijk in +038 Ki, 


berekend met een 


km (voor een 


tuurpiek in het 


3 ba an nt En N ie . EA zen - N « za En 
Later overschakeling op rooster 39 x 39 met DX = 2019 km en DZ = 42/5 Km 


hd 


op Tr = 1252 sec. De nieuwe tussenliggende punten gedefinie 


Begintoestand: drukevenwicht 


lop i5 x 15 rooster) 


rd als 





centrale 9 punten : dichtheid m j.l x 10° em ® 

dearomheen: s id « 1,08 & 10° em ? 
ä ; E = 1.04 «10° em ° 
k : 8 1.0 & 10° cm ° 
k 8 1.003 «10° em * 
k 4 = 1.001 # 10° em ° 
ì ; 0 = 1.000 # 10° cm ° 


Mechanische verhitting *“ ps aan rand geen nettovenergieverlies, 


a. dichtheidscontcuren op t = 1252 sec (1 « dichtheid 1.01 « 10° cm %; 


Ve B 


10 = dichtheid 1.47 x 10° em °). 
b, dichtheidscontouren op t = 1580 sec (1 = dichtheid 1,06 « 10° cm °5 
: IO = dichtheid 3.83 « 10° cm °). 


io 
a 


dichtheidscontouren op t = 1760 sec (ì « dichtheid 0.94 & 10° cr °3 
10 = dichtheid 22,5 « 10° cm °) | 
snelheidsvektorplot op t = 1580 sec; maximale snelheid = 11.5 km/serc; 


f 
: Ga 
Ee 


de kleinste snelheidsvektoren duiden aan een 100 keer kleinere 
snelheid aan dan de grootste. 


2e snelheidsvektorplot op t = 1760 secs maximale snelheid = 209 km/sec. 


Figuur 3: 
Finch effekt, 
| : | Tweedimenstionaal roosters 2i x 15 punten. 
| x=richting (verticaal); DX = 7778 km. 
zerichtiag (horizontaal); DZ = 2333 km. 
Cox & Tucker stralingsfunctie, 
Geen warmtegeleiding. 
Begincendities: gebied met 10% dichtheidsperturbatie met afmetingen 
(7 DX) x (13 DZ) = 5500 x 30 000 km. | 
Begindichtheid buiten geperturbeerd gebied:p =5.8RIO 1fo/cn. 
Te tf K 
Mechanische verhitting "“ dichtheid; aan rand geen netto= 
energieverlies. 
a. dichtheidsprofielen op verschillende tijdstippen: 
{0) Cm OQ sec 
(1) t = 1300 sec. 





od 


ET. ES 


(2) t = 2399 zee 
(3) t m= 2700 sec 
(4) t = 2800 sec 
(53 £ « 3900 sec 
(6) t  m 2953 sec 
b. Breedte (in de z-richting) van het kaale, dichte gebied als fvactie 
van de tijd uitgedrukt in de staplengte D4, 
&. Temperatuur als functie wan de tijd voor verschillende punten langs 
de z=ag8, 
O + 1 = punt in centrum + punt naast centrum (op de z-as). 
è m twee posities vanuit het centrum langs de 2786, 
UZ o 
d, snelbeidsprofielen (langs de z-as) op verschillende tijdstippen 
(1) = 2996 sec. 
(8) = 3046 sen. 
Thermische instabiliteit, met en zonder warmtegeleiding. 
Cox & Tucker functi | 
Zl x 15 punten, 
DX = 7773 km; DZ » 2353 km. 
Mechanische verhitting evenredig met dichtheid; aan rand geen netto 


energieverlies, 


Begincondities: gebied wet 10% dichtheidsperturbatie met afmetingen 


(3DX) x {3DZ) "* 22000 x 7000 km. 


Begindichtheid buiten geperturbeerd gebied p‚=5.Ba10 TPg/en’ 
Tm 10° KR 
e 


De druk is overal Ì. 


netto energieverlies per cm° per sec (RAD) en temperatuur als functie 


‚ 


van de tijd. 


TEMP, = temperatuur als functie van de tijd zonder warmte 


gelel idiug. 


TEMP, = temperatuur als functie van de tijd met een isotrope warmte” 


geleidingscoefficiënt 0.01 « K 


5 
ĳ 


corona zonder B-veld. 


nette energieverlies (per sec per cm°) zonder warmtegeleiding. 
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hes 


RAD, » idem, nu met warmtegeleiding. 

Db. deuk (p) en dichtheid (RHO) als functie van de tijd met en zonder 
warmtegeleiding. | 
RHO ja per definitie 1,1 op t = Q in het geperturbeerde gebied, de 
druk is dan Ì. 

C. enkele teuperatuurprofielen op tijdstippen waarop de invloed wan de 
warmtegeleiding van dezelfde orde van grootte is ale de stralings= 
afkoeling. OO | nn 

d. temperatuurprofielen op tijdsti ppen waarop stralingsafkoeling over= 
heerst, 

e. _RÀD=-profielen. 


£. Arukprofielen. 


Ziguur Hi: 

Dezeifde gegevens als figuur 10, zonder warmtegeleiding, alleen nu berekend 
met de McWhirter stralingsfunctie. 

a. temperatuur als functie van de tijd en RAD als functie van de tijd. 

hb. dichtheid (RHO) en druk (p) als functie van de tijd. | 


Figuur 12: 

Ket netto energieverlies per sec als functie van de temperatuur voor de 
Cox & Tucker eu de McWhirter functie voor een 10% perturbatie voor de 
dichtheid, | 

8. voor een Zeoberisch proces. 

b. zoals berekend met het 21 x 15 rooster van £ig.î0 en 11. 

RAD is hier Ì voor de Cox & Tucker functie als de temperatuur 10°R/S.3 
bedraagt. Dan is RAD woor de McWhirterfunctie G.i85, 


Ke 


De mechaaische verhitting ie evenredig met de dichtheid verondersteld. 


Figuur 13: 


al | 
x wassaflux uit: t-x=richtin , en 
De verhouding == = ngen als functie van de tijd, woor 
M_ massafiux uit t z-richting 


de gevallen zonder warmtegeleiding van figuren 10 en Ii (zie vok figuur 


cp blz.15 en tekst bla.21). 
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